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Inorganic perovskite oxide ferroelectrics have recently generated substantial interest for photovoltaic 
applications; however, existing materials suffer from excessive electronic band gaps and insufficient 
electric polarization. The recent resurgence of hybrid perovskite ferroelectrics addresses the afore¬ 
mentioned deficiencies, but they are highly unstable against environmental effects and an inorganic 
resolution may still be optimal. Here we use first-principles calculations to design a low band gap, 
planar ferroelectric by leveraging the complexity of layered double perovskite oxides AA'BB'Oe- We 
exploit AylA' size mismatch, a nominally empty d-shell on the B-site, and A cations bearing lone- 
pair electrons to achieve a large ferroelectric polarization. Additionally, B' ^B is chosen to achieve 
full charge transfer and a Mott susceptible filling on the B'-site, yielding a low band gap. These 
principles are illustrated in BaBiCuVC>6, BaBiNiVOe, BaLaCuVOe, and PbLaCuVOe, which could 
be realized in layer-by-layer growth. This new class of materials could lead to stable, high-efficiency 
photovoltaic. 


I. INTRODUCTION 

Ferroelectric oxides have emerged as a candidate for 
applications in photovoltaic devices, which separate 
photo-excited electron-hole pairs via the internal electric 
field of the ferroelectric 1 . Recent progress in domain en¬ 
gineering of ferroelectric oxides has demonstrated above- 
bandgap open circuit voltages 2,3 , in contrast to classic 
p-n junction photovoltaic 4 . Furthermore, ferroelectric 
oxides are tolerant to various environments and can be 
fabricated using cost efficient methods 5,6 . However, their 
power conversion efficiencies (PCE) have been quite lim¬ 
ited until recently, despite decades of exploration in this 
class of materials 5,7 9 . The limitation of many existing 
ferroelectric oxides is a large electronic band gap (ie. > 
3eV) 7 , which is typical of do systems where the gap is be¬ 
tween the oxygen 2 p and the nominally empty transition- 
metal d states. Band gaps of this magnitude are prob¬ 
lematic, as they cannot capture a significant component 
of the visible spectrum 4 . An additional limitation is the 
extremely low photo-current observed in most ferroelec¬ 
tric oxides 7 , and this problem was long hindered by a 
lack of a detailed mechanistic understanding. However, 
the application of shift current theory has provided sig¬ 
nificant insight into this phenomenon, and it is apparent 
that strong polarization and covalency are necessary in¬ 
gredients for enhancing the photocurrent 10,11 . 

In order to circumvent the aforementioned limitations, 
one must simultaneously engineer a low band gap, a large 
ferroelectric polarization, and strong hybridization. In¬ 
deed, designing a lower band gap in ferroelectric oxides 
has drawn a large degree of attention 5,8,12-14 . One possi¬ 
ble route to smaller band gaps is to utilize Mott or charge- 
transfer insulators, as in the well-known case of BiFeOs 12 
where an optical band gap of 2.7eV is realized. While 
the gap is still excessive in this case, this is a promis¬ 


ing direction that can be furthered in the context of al¬ 
loys. In particular, one can explore the double perovskite 
A 2 BB'0 6 in order to combine both aspects of band insu¬ 
lators (ie. nominal do valence) with Mott insulators (ie. 
nominal d n valence) to achieve a high degree of tunability 
of the band gap. A large ferroelectric polarization in the 
presence of strong hybridization can be achieved via the 
second order Jahn-Teller effect, as is often realized when 
A-site cations contain lone-pair electrons and or when B- 
site ions are in a nominal do valence. Furthermore, the 
magnitude of the ferroelectric polarization can be fur¬ 
ther enhanced by employing size mismatch on the A/A'- 
sites, particularly when one of the cations is lone-pair 
active 15-18 . In order to maximally exploit all aforemen¬ 
tioned mechanisms, one can consider the more general 
double perovskite AA'BB 7 06 19-21 . Such compounds can 
be realized in experiment via chemical synthesis 5,22 24 
and potentially via layer-by-layer growth such as Molec¬ 
ular Beam Epitaxy (MBE) and Pulsed Laser Deposition 
(PLD) 9,25,26 . One prominent experiment was the chemi¬ 
cal synthesis of [KNbOsJi-^BaNq^Nb^Os-^, with 
x ~ 0.1 , whereby the large band gap of KNbOs (ie. 
3.8 eV) is reduced to 1.39 eV; though the polariza¬ 
tion is reduced from 0.55Cm ~ 2 to 0.19 Cm ~ 2 , the pho¬ 
tocurrent is still ~ 50 times that of the classic ferroelec¬ 
tric (Pb,La)(Zr,Ti)03 5 . Another example with A=A' is 
Bi 2 FeCr06, resulting in a band gap of 1.43 - 2.51 eV, 
depending on cation ordering, and a ~8Yo power con¬ 
version efficiency (PCE) 9 . However, there is also a sig¬ 
nificant reduction of ferroelectric polarization, which is 
detrimental for photovoltaic performance. The preced¬ 
ing examples illustrate the potential of enhancing photo¬ 
voltaic performance via band gap reduction, though the 
ferroelectric polarization was excessively suppressed in 
these cases, which likely degrades the PCE. Therefore, 
first-principles approaches can be a powerful avenue to 
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more broadly explore this vast phase space, and design 
low band gap materials which do retain large ferroelectric 
polarization with strong hybridization. 

Engineering double perovskites from first-principles 
has received some attention within the literature. Previ¬ 
ous studies have explored BB' combinations resulting in 
do/dio valence states, but these scenarios always result 
in large band gaps 16,17 ’ 27 . Expanding to Mott physics, 
d n /dm scenarios have been successful in attaining low 
band gaps 15 but experimental results reveal relatively 
small polarizations (ie. 0.2 Cm~ 2 ) and the origin of 
ferroelectricity is AA' size mismatch which is not as fa¬ 
vorable in terms of hybridization 28 . Another promising 
direction is BB' combinations resulting in do/d n valen¬ 
cies, which may result in low band gaps 29 . However, 
it has been suggested that the nature of the BB' or¬ 
dering will be critical to realizing a polarization asso¬ 
ciated with the do ion 24 , with rock salt ordering being 
highly unfavorable. This argument is consistent with 
the lack of polarization in rock salt ordered double per- 
ovskite A 2 MWO 6 (M=Ni,Mn,Cu,etc) as determined in 
experiment 30-33 .Therefore, it is desirable to pursue lay¬ 
ered ordering of BB'. Furthermore, octahedral tilting 
competes with the second order Jahn-Teller distortion as¬ 
sociated with the do site 34 , and having a layered ordering 
of AA' with large size mismatch can frustrate octahedral 
tilting 35 . 

In this letter, we design a new class of materials where 
we utilize all of the aforementioned mechanisms simulta¬ 
neously to achieve a low electronic band gap and a large 
ferroelectric polarization with strong hybridization. We 
employ at least one A-type ion with lone pair electrons 
(ie. Bi 3+ or Pb 2+ ). In order to maximize AA' size mis¬ 
match, the A' ion is chosen to have a different valence 
than the A-type. In particular, Bi 3+ is paired with Ba 2+ 
while Pb 2+ is paired with La 3+ , where the selection was 
dictated by maximizing the size of the A-type ion in the 
respective valence state. We realize full charge transfer 
with select BB' pairs having valence d n /do, which results 
in a low band gap. Given our A A' selections, nominal 
valence counting dictates (BB') +7 , and one of these ions 
should be do- Therefore, we select vanadium as the do 
ion, which must exist in the 5 + valence, in conjunction 
with various B ions in the 2 + valence. Vanadium is cho¬ 
sen over, for example, Ti 4+ or Nb 5+ for several reasons. 
First, V 5+ is smaller than Nb 5+ and Ti 4+ , 36 which will 
lead to a larger polarization. Second, having a do B' 5+ ion 
dictates a nominal B 2+ -type ion that is more susceptible 
to being a Mott insulator for the later transition metals, 
which are necessary in order to ensure full charge trans¬ 
fer (ie. sufficient electronegativity differences). There¬ 
fore, for B-type ions we choose Ni 2+ and Cu 2+ . These 
physically motivated choices are permuted to search for 
optimal properties, and we additionally explore the va¬ 
lidity of these underlying rules. A final condition is that 
we restrict this study to layered ordering of AA' and BB' 
ions, due to the potential for realization via layer-by- 
layer growth, such as MBE and PLD. We begin with a 


detailed exploration of BaBiCuVOe, which will serve as 
our prototype, and then generalize to other possible per¬ 
mutations. 


II. METHOD 

We utilized density functional theory calculations im¬ 
plemented in the Vienna ab initio simulation program 
(VASP) 37-40 , which employs a plane-wave basis and 
projector-augmented wave (PAW) potentials 41,42 . Cor¬ 
relation effects on the 3d orbitals are treated within the 
local spin density approximation (LSDA) + U method 43 . 
The LSDA + U method requires the definition of the on 
site Coulomb repulsion V, and we do not include an on 
site exchange interaction J given that this is accounted 
for within LSDA 44 . We explore a wide of range of values 
for I/, given that this is still an open methodological prob¬ 
lem. All results were obtained using 6x6x6 Monkhorst 
- Pack k-point mesh centered at T and a 500 eV plane 
wave cutoff 45 . The tetrahedron method is used for Bril- 
louin zone integration 46 . In order to find the ground 
state structure, we have explored several octahedral ro¬ 
tation patterns and computed phonon dispersion curves. 
Phonons were computed using the PHONOPY code and 
a 2 x 2 x 1 supercell with forces obtained from self- 
consistent total energy calculations using the finite dis¬ 
placement method. 47 The lattice parameters and atomic 
positions were relaxed at the LSDA+U levels until the to¬ 
tal energy changed by less than 10~ h eV per unit cell and 
the residual force were smaller than 0.01 eV per A. Born 
effective charges and the ferroelectric polarization were 
computed within the modern theory of polarization 48,49 . 
The endpoints of the ferroelectric path are defined us¬ 
ing a 180° in-plane rotation centered at the vanadium 
atom. Nudged elastic band is performed to determine 
energy barrier of the ferroelectric path 50 . Ionic Radii are 
given by Shannon for different coordination and charge 
states 36 . All images of crystal structures were generated 
using VESTA code 51 . 

III. RESULTS 

A. Reference structures and spontaneously broken 
symmetries 

In order to elucidate the origin and nature of the fer¬ 
roelectric polarization in layered AA'BB'06, we consider 
three reference structures. The first reference structure, 
referred to as R 0 , consists of placing the atoms on the 
ideal perovskite lattice and then allowing full relaxations 
of the lattice parameters while only relaxing the atoms in 
the z-direction (space group Cmm2, see Figure la). The 
second reference structure, R q =ce corresponds to fully 
relaxing R 0 with respect to lattice parameters and all in¬ 
ternal coordinates within the unit cell (space group Cm , 
not pictured). The third reference structure, 
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corresponds to creating a supercell with lattice vectors 

(1/V2,1/V2,0),(1/V2,-1/V2,0) and (0,0,1) in 

terms of R 0 , and then fully relaxing lattice parameters 
and all internal coordinates (space group Pc, see Fig¬ 
ure lb). These structures allow one to determine the 
role of q = 0 relaxations and octahedral rotations. It 
should be noted that even the highest symmetry refer¬ 
ence R 0 contains a static polarization in the z-direction, 
though an in-plane ferroelectric polarization can still be 
formed by spontaneously breaking C± v symmetry via the 
destruction of a mirror plane upon degenerating to P q= o 

Or Pq=7r,7r,0* 

We begin by examining the R 0 structure for 
BaBiCuVOe (see Figure la, lb). In this case, we observe 
a large structural relaxation of both the Ba and Bi to¬ 
wards the CU-O 2 layer, presumably due to the large nom¬ 
inal 5+ charge associated with the V ion. This strong dis¬ 
tortion will be seen in all subsequent structures, irrespec¬ 
tive of octahedral rotations and other further symmetry 
breaking. While there is a polarization in the z-direction 
associated with this A/A' ordering and the subsequent 
z-direction distortion, this is not relevant for our inter¬ 
ests as this ordering is static and could not practically be 
switched. The main question is therefore how symmetry 
is spontaneously broken relative to this reference. 

In figure 2, we have computed phonon spectra for the 
three reference structures of BaBiCuVOe (see Methods 
Section). In the R 0 structure there are many soft-modes, 
including at the T-point, reflecting the unstable nature 
of this highest symmetry structure (see Figure 2a). In 
the R q =o structure, all of the T-point instabilities have 
been removed as this structure is a fully relaxed degen¬ 
eration of the R 0 structure (see Figure 2b). This effect 
is pictorially illustrated in the comparison of the planar 
projections of Figures 2f and 2g. However, there are still 
strong in-plane instabilities shown in the phonons, which 
correctly suggests that there is some lower energy tilting 
pattern of the octahedron. In the R q=7r ^ 0 structure, it 
is demonstrated that all in-plane instabilities have been 
removed (see Figure 2c), and the resulting structure is 
illustrated in the planar projections of Figure 2h and 2i. 
However, the R q=7r ^ 0 structure still has a small instabil¬ 
ity in the z-direction of reciprocal space. By doubling the 
unit cell in the z-direction and allowing for full atomic re¬ 
laxations, we have determined that this instability only 
causes a very small structural distortion that does not 
substantially affect the physics (ie. the energy gain is 
~ 1 meV per f.u.). Therefore, R q=7r ^ 0 is a sufficient ap¬ 
proximation for our ground state structure. 

Having established how the symmetry degenerates, we 
now use the Born effective charges to give an approxi¬ 
mate layer-resolved ferroelectric polarization, referenced 
by R 0 . The main contribution to the polarization comes 
from large displacements within the Bi-0 and V-O 2 lay¬ 
ers, which arise from the second order Jahn-Teller effect 
associated with the lone-pair electrons and do configu¬ 
ration, respectively. In the R q= o structure, the polar¬ 
ization is P = 0.435i + 0.237j Cm -2 , demonstrating 
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FIG. 1 . | Reference structures for layered BaBiCuVOe. 

(a) Reference structure R 0 , where the atoms have been 
placed on an ideal perovskite lattice and subsequently re 
laxed only in the z-direction. As illustrated, this structure 
retains symmetry, though it is unstable, (b) Refer¬ 
ence structure R q = 7 r, 7 r,o, which allows for ful 1 relaxation of 
all q =s (7r, 7r, 0) distortions and lattice strains. 


that a large polarization is formed even in the absence 
of any possible octahedral tilts. The layer-resolved anal¬ 
ysis further provides the polarization of Bi-0 layer and 
V-O 2 layer, resulting in P = 0.248i + O.llOj Cm -2 
and P = 0.165i + O.llOj Cm 2 , respectively. In the 
R q=7r?7r 0 structure, the magnitude of the polarization 
remains similar though the direction changes yielding 
P = 0.354i + 0.354j Cm -2 , which can further be split 
into P = 0.189i + 0.189j Cm -2 for the Bi-0 layer and 
P = 0.127i + 0.127j Cm -2 for the V-O 2 layer. There¬ 
fore, octahedral tilts do perturb the direction of the po¬ 
larization. It should be noted that the polarization in 
the Bi-0 layer is aligned in the same direction as the po¬ 
larization in the V-O 2 layer, giving rise to a very large 
in-plane polarization. The possibility of antialigning the 
Bi-0 and V-O 2 layers was explored, but this could not 
be stabilized in our calculations. 


B. Energy barriers and polarization 

In order to achieve a precise understanding of the fer- 
roelectricity, we employed the modern theory of polar¬ 
ization to study the ferroelectric moment of BaBiCuVOe 
in more detail. The modern theory of polarization has 
theoretically reproduced the precise magnitude of the fer¬ 
roelectric moment in many systems 48,49 . In our system, 
spontaneous breaking of C± v symmetry inherently pro¬ 
vides disparate paths between symmetry equivalent min¬ 
ima, and we explore the paths defined by the Rq =7r7r o 
structure rotated by C^ and C^/ 2 , respectively. To be¬ 
gin with, the path defined by C^ is computed using the 
modern theory of polarization yielding 0.65 Cm~ 2 (see 
Figure 3a). This value is in the vicinity of the polar- 
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FIG. 2. | Phonons of successive reference structures 

The phonons of reference structures for BaBiCuVOe are pre¬ 
sented in panel (a) ,(b), (c) corresponding to R 0 , R q =o, and 
Rq= 7 r, 7 r,o, respectively. Correspond ing in-plane projections of 
the Bi-0 and V-0 layers are given in panels (d) ,(e) for R 0 , 
panels (f), (g) for R q =o, and panels (h), (i) for R q= 7 r,7r,o- 
In successive order, the polarization is 0, 0.435i + 0.237j 
Cm ~ 2 , and 0.354i + 0.354j Cm~ 2 . 


ization of BiFeOs 52 , which is 0.8 Cm~ 2 . Additionally, 
the corresponding energy barrier is calculated to be 277 
meV per transition metal ion (see Figure 3b). Alter¬ 
natively, the path defined by C ^/2 results in half of the 
energy barrier, and this could lead to enhanced switching 
speeds 53 . These energy barriers for BaBiCuVOe are rea¬ 
sonable in the sense that BiFeOs has a theoretical energy 
barrier of 427 meVperf.u. (ie. within Generalized gradi¬ 
ent approximation + U) yet the ferroelectric moment is 
experimentally switchable 54 . 

The differences in total polarization between the 
Berry-phase approach and Born effective charge analysis 
reflects the presence of anomalous dynamical charges 55 
in BaBiCuVOe- In other words, it highlights the cova¬ 
lent character of the atomic bonds and, in turn, of the 
electronic structure of BaBiCuVOe, which is also seen 
in BaTiOs 55 . Additionally, the born effective charges of 
active ions in the R 0 structure exhibit a striking devi¬ 
ation from nominal charges, which clearly demonstrates 
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FIG. 3. Energy barrier and ferroelectric moment | (a) 

Total polarization P, (b) total energy profile E tot as a func¬ 
tion of the polar distortion characterized by A (ie. A = 0 
denotes the ferroelectric phase and A = 1 is the midpoint). 
The energy barrier is obtained via the nudged elastic band 
method and polarization is calculated through both the mod¬ 
ern theory of polarization and Born effective charge analysis 
of fully distorted structure 

a highly covalent character. More specifically, using the 
Born effective charges from the R 0 structure yields a po¬ 
larization of 0.81 Cm ~ 2 , while using the Born effective 
charges from the R q=7r ^ 0 structure yields 0.54 Cm~ 2 . 
The average of those two values is similar to the value 
from the modern theory of polarization, 0.65 Cm~ 2 . 
These observations are a clear demonstration of the im¬ 
portance of hybridization effects. 


C. Electronic structure and magnitude of band gap 

Here we explore the basic electronic structure of 
BaBiCuVOe, which follows the simple rules outlined in 
the introduction. More specifically, V nominally donates 
its sole d-electron to Cu, analogous to what was observed 
for Ti in La2NiTi06 29 . We begin by illustrating this via 
density functional theory (DFT)+U calculations, where 
we use a Hubbard U of 5 eFon the copper site, consis¬ 
tent with previous studies 56 , and we conservatively use 
U =0 on the V site. The effect of the on-site Hubbard U 
will be explored in detail below. The standard fully lo¬ 
calized limit double counting 57,58 is used, and no on-site 
exchange is employed given that we are using a spin- 
dependent density functional 44 . The resulting density of 
states (DOS) for the ground state antiferromagnetic spin 
ordering yields a band gap of 1.1 eV (see Figure 4a). The 
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FIG. 4. | Electronic structure Atomic resolved density of 
states of BaBiCuVOe in the Rq =7 r, 7 r,o structure for (a) Uv 
— 0 and Ucu — 5 ; (b) Uv = 0 and Ucu = 0 eV ; (c) 
projection onto copper d-orbitals; Fermi energy is set to be 0; 
(d) The variance of band gap depending on Uv and Ucu- 


atom resolved DOS for V exhibits a strong peak above 
the gap indicating a complete nominal charge transfer, 
which is corroborated by zero magnetic moment on the 
V atom. However, there is a substantial orbital occupa¬ 
tion of ~ 4 electrons for the V d-orbitals due to strong 
hybridization with the oxygen, as illustrated by the sub¬ 
stantial DOS projection of V towards the bottom of the 
oxygen bands. Below we demonstrate that this results in 
the band gap being very insensitive to the value of U on 
the V site. Cu has a magnetic moment of 0.5 ub, which 
is substantially reduced from the nominal value of 1 he- 
Projecting the DOS onto the Cu-3d states reveals states 
above the gap, which are driven by the on-site Coulomb 
repulsion U (see Figure 4a, 4b). This scenario is analo¬ 
gous to the cuprates, where there is a substantial crystal 
field that breaks E g symmetry and results in an effective 
single band Hubbard model which yields an antiferro¬ 
magnetic Mott insulator. This is clearly illustrated by 
the orbitally resolved DOS presented in Figure 4c, where 
d xy character is seen both above and below the gap while 
d z 2 is nearly completely filled. 

The precise value of the on-site Coulomb repulsion U is 
still an open problem, and the DFT+U method is crude 
in the sense that it is a Hartree approximation to the 
dynamical mean field theory(DMFT) impurity problem 
within the DFT+DMFT formalism. Both of these issues 
warrant a careful exploration of how the results are af¬ 
fected by changes in U. Therefore, we compute the band 
gap of the system as a function of U for vanadium and 
copper, exploring essentially all feasible values (ie. 0 - 
6eV, 0 - 10 eV for vanadium and copper, respectively; 
See Figure 4d). The U of vanadium has very little effect 
due to the large hybridization with the oxygen. Neglect¬ 
ing exchange, the orbital dependent term within DFT-fU 
that is added to the Kohn-Sham potential has the form 
U( 1/2 — n*) 43 , such that orbitals near half filling will not 
be substantially perturbed unless U is sufficiently large 
relative to the hybridization. Alternatively, the U on Cu 
is essential for driving the Mott physics and is necessary 
to open a gap (see Figure 4a, 4b). Increasing the value 
of V increases the gap, as expected, though the magni¬ 
tude of the increase is relatively small due to the strong 
hybridization of Cu and oxygen, and eventually the band 
gap largely saturates due to the upper band edge being 
dominated by V. Even for a relatively large value of U = 
6 eV, the band gap remains relatively small at 1.35 eV. 
This scenario should be contrasted to La2TiNi06, where 
the gap depended strongly on the U value for Ti 29 due to 
a relatively smaller electronegativity difference for Ti-Ni 
and less Ti-0 hybridization. 


D. Generalization to other compounds 

While we have illustrated our design principles in the 
context of BaBiCuVOe, these rules can obviously be 
used to generate numerous other compounds with sim¬ 
ilar properties. To illustrate this point, we have sue- 
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TABLE I. The summary of planar Ferroelectric Mott | The electronic gap, energy barrier, (layer-resolved) polarization, 
and AA' size difference based on Shannon radii 36 for AA'BB'Oe compounds. The energy barrier is normalized per transition 
metal. Total polarization is computed by modern theory of polarization. The asterisk indicates that U was increased to 6.5 
eV and 8.0 eV on Ni and Cu, respectively, in order to ensure a gap was maintained over the entire distortion path, which is a 
necessary condition for the modern theory of polarization. This change has a negligible effect on the other observables. 


AA'BB'Oe 

BaBiCuVOe 

PbLaCuVOe 

BaBiNiVOe 

BaLaCuVOe 

Gap (eV) 

Energy barrier (meV) 
Polarization (Cm -2 ) 
A/AAize mismatch ( A ) 
Unit Cell Volume (A 3 ) 
A-O (Cm -2 ) 

A'-O (Cm -2 ) 

B-O (Cm -2 ) 

B'-O (Cm -2 ) 

1.08 

287 

0.65 

0.25 

240.47 

0.04 

0.28 

0.00 

0.19 

1.00 

226 

0 . 52 * 

0.13 

229.72 

0.12 

0.10 

0.00 

0.23 

1.05 

269 

0 . 68 * 

0.25 

237.41 

0.03 

0.28 

0.03 

0.22 

0.86 

214 

0 . 41 * 

0.26 

239.38 

0.03 

0.08 

0.00 

0.16 


cessfully realized three additional candidates in addi¬ 
tion to BaBiCuVOe: PbLaCuVOe, BaBiNiVOe, and 
BaLaCuVOe (see Table 1). We first consider the case 
of PbLaCuVOe, where we have decreased the A/A' mis¬ 
match while retaining one lone-pair cation and kept ev¬ 
erything else identical. As expected 16 , the smaller A/A' 
mismatch results in a smaller polarization within the 
lone-pair cation layer, though a nontrivial polarization 
does form in the La-0 layer. Nonetheless, the net ef¬ 
fect is an overall decrease in polarization and a large 
decrease in the energy barrier, which could be benefi¬ 
cial for ease of electrical switching. Due to the relatively 
small changes in the CU-O 2 planes, the band gap changes 
very little. Another interesting possibility is to retain 
the Ba/Bi cations and instead exchange Ni in place of 
Cu. In this case, Ni will result in a 2+ state, which is 
susceptible to Mott physics, similar to NiO and the sce¬ 
nario outlined in La2TiNi06 29 . We find a similar barrier 
and polarization to BaBiCuVOe, which is consistent with 
the only change being associated with the ion exhibiting 
Mott physics. Additionally, the band gaps are relatively 
similar, assuming a reasonable value of U =5eV for Ni 59 . 
The final scenario we explore corresponds to the removal 
of lone-pair cations, which is the case of BaLaCuVOe- 
Here the only nontrivial polarization in the A/A' layers 
is associated with La-O, which is analogous to the large 
size mismatch as in the KNbOs-LiNbOs system 60 . The 
band gap is slightly reduced relative to BaBiCuVOe due 
to diminished tilting in the CU-O 2 layer. 


IV. DISCUSSION 

Our proposed class of planar ferroelectric Mott insu¬ 
lators optimally matches the criteria for use as a photo- 
ferroelectric 4,61 : band gap of ~ leV, strong ferroelectric 
polarization and covalent character, and a large phase 
space of chemical elements such that properties can be 
tuned. The low band gap is essential to absorb a sig¬ 
nificant portion of the solar spectrum. The sizable po¬ 
larization and strong covalent character are necessary to 
obtain a substantial shift current 10,11 , which is believed 
to be indicative of the photocurrent and hence the power 


conversion efficiency (PCE). Our current work has iden¬ 
tified four compounds which satisfy our design rules, and 
there are many more compounds which could be included 
to form a substantial phase space of possibilities. Addi¬ 
tionally, the fact that our systems are antiferromagnet 
Mott insulators could further enhance the potential of 
these materials in two ways 62 . First, multiple charge car¬ 
riers through impact ionization suggest the possibility of 
exceeding the Shockley-Queisser efficiency limit 4 . Sec¬ 
ond, high drift voltage and large diffusion length could 
be helpful despite a low mobility due to strong correla¬ 
tions. Antiferromagnetic Mott insulators have demon¬ 
strated the possibility of enhancing ultrafast separation 
of photo doped carriers before recombination 63 . Only 
careful experiments will be able to definitively determine 
the applicability of Mott insulators in this context due to 
the complexity of strongly correlated electron systems. 

Another potential area of application exists in the con¬ 
text of non-volatile memory. The photovoltaic output of 
BiFeOs was employed as a read-out signal, switching be¬ 
tween 0.11 V and -0.23 V for up to 108 cycles 64 . In 
terms of non-volatile memory devices, an in-plane fer¬ 
roelectric polarization has several benefits 2 . First, in¬ 
plane ferroelectrics should be much less sensitive to ex¬ 
plicitly broken symmetry induced by a substrate, lead¬ 
ing to a more symmetric voltage and current bi-state. 
Second, in-plane ferroelectrics may be straightforwardly 
exploited via the growth of large area thin-films gener¬ 
ating a large voltage, which would be difficult to achieve 
in a thick film out-of-plane ferroelectric. An open cir¬ 
cuit voltage of 8 volts is achieved with a small inci¬ 
dent intensity of 1 mWcm~ 2 for a thin film device of 
Pbo.97Lao.o3(Zro.52Tio.4803. Thus our proposed class of 
materials could be multifunctional in the sense that the 
ferroelectric order is coupled to optical and electronic 
transport channels. This could lead to optically ad¬ 
dressed ferroelectric memory with non-destructive read 
out 65 . It is also worth noting that a large planar ferro¬ 
electric moment could assist the formation of charged 
domain walls, which have their own potential device 
applications 66 . 

It is worthwhile to compare our systems with the re¬ 
cent resurgence of organometallic halides in the context 
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of photovoltaic, which has also resolved the aforemen¬ 
tioned deficiencies of ferroelectric oxide photovoltaic 61 . 
In the organometallic halides, variation in hysteresis of 
the current-voltage profile under different scanning con¬ 
ditions implies ferroelectric behavior and exceptionally 
long carrier diffusion lengths, leading to the hypothesis 
of “ferroelectric highways”. Recent ab-initio calculations 
confirm a large polarization and strong covariance in pho- 
tovoltaicly promising hybrid lead halide perovskites 67,68 , 
which is consistent with the general view of shift current 
theory 10 . Nonetheless, their instability against water due 
to weak halide bonding and environmental issues associ¬ 
ated with lead are still very serious problems 69 , and re¬ 
solving these deficiencies could be formidable. Therefore, 
the approach outlined in this paper could be a very valu¬ 
able path to explore in parallel with the recent trends in 


organomet allies. 


ACKNOWLEDGMENTS 

We would like to thank D. Vanderbilt and K. Rabe for 
valuable discussions. The authors acknowledge funding 
from the FAME, one of six centers of STARnet, a Semi¬ 
conductor Research Corporation program sponsored by 
MARCO and DARPA. The research used resources of the 
National Energy Research Scientific Computing Center, 
a DOE Office of Science User Facility supported by the 
Office of Science of the U. S. Department of Energy under 
Contract No. DE-AC02-05CH11231 


1 V. M. Fridkin, Photoferroelectrics (Springer-Verlag, 
Berlin, 1979). 

2 K. Yao, B. K. Gan, M. Chen, and S. Shannigrahi, Applied 
Physics Letters 87 , 212906 (2005). 

3 S. Y. Yang, J. Seidel, S. J. Byrnes, P. Shafer, C.-H. Yang, 
M. D. Rossell, P. Yu, Y.-H. Chu, J. F. Scott, J. W. Ager, 
L. W. Martin, and R. Ramesh, Nature nanotechnology 5, 
143 (2010). 

4 W. Shockley and H. J. Queisser, Journal of Applied Physics 
32, 510 (1961). 

5 I. Grinberg, D. V. West, M. Torres, G. Gou, D. M. Stein, 

L. Wu, G. Chen, E. M. Gallo, A. R. Akbashev, P. K. 
Davies, J. E. Spanier, and A. M. Rappe, Nature 503 , 
509 (2013). 

6 D. Cao, C. Wang, F. Zheng, W. Dong, L. Fang, and 

M. Shen, Nano letters 12 , 2803 (2012). 

7 Y. Yuan, Z. Xiao, B. Yang, and J. Huang, J. Mater. Chem. 
A 2 , 6027 (2014). 

8 G. Zhang, H. Wu, G. Li, Q. Huang, C. Yang, F. Huang, 
F. Liao, and J. Lin, Scientific reports 3, 1265 (2013). 

9 R. Nechache, C. Harnagea, S. Li, L. Cardenas, W. Huang, 
J. Chakrabartty, and F. Rosei, Nature Photonics 9, 61 
(2014). 

10 S. M. Young and A. M. Rappe, Physical Review Letters 
109, 116601 (2012). 

11 S. M. Young, F. Zheng, and A. M. Rappe, Physical Review 
Letters 109, 236601 (2012). 

12 T. Choi, S. Lee, Y. J. Choi, V. Kiryukhin, and S.-W. 
Cheong, Science (New York, N.Y.) 324 , 63 (2009). 

13 W. S. Choi, M. F. Chisholm, D. J. Singh, T. Choi, G. E. 
Jellison, and H. N. Lee, Nature communications 3, 689 
( 2012 ). 

14 S. R. Basu, L. W. Martin, Y. H. Chu, M. Gajek, 
R. Ramesh, R. C. Rai, X. Xu, and J. L. Musfeldt, Ap¬ 
plied Physics Letters 92 , 091905 (2008). 

15 D. J. Singh and C. H. Park, Physical Review Letters 100, 
087601 (2008). 

16 S. Takagi, A. Subedi, D. J. Singh, and V. R. Cooper, 
Physical Review B 81 , 134106 (2010). 

17 S. Takagi, A. Subedi, V. R. Cooper, and D. J. Singh, 
Physical Review B 82 , 134108 (2010). 

18 T. Shishidou, N. Mikamo, Y. Uratani, F. Ishii, and 


T. Oguchi, Journal of Physics: Condensed Matter 16, 
S5677 (2004). 

19 J. Young, A. Stroppa, S. Picozzi, and J. M. Rondinelli, 

Dalton transactions (Cambridge, England : 2003) 44, 

10644 (2015). 

20 J. M. Rondinelli and C. J. Fennie, Advanced materials 
(Deerfield Beach, Fla.) 24, 1961 (2012). 

21 A. T. Mulder, N. A. Benedek, J. M. Rondinelli, and C. J. 
Fennie, Advanced Functional Materials , n/a (2013). 

22 P. Zubko, S. Gariglio, M. Gabay, P. Ghosez, and J.-M. 
Triscone, Annual Review of Condensed Matter Physics 
( 2011 ). 

23 T. Fukushima, A. Stroppa, S. Picozzi, and J. M. Perez- 
Mat o, Physical chemistry chemical physics : PCCP 13, 
12186 (2011). 

24 M. C. Knapp and P. M. Woodward, Journal of Solid State 
Chemistry 179, 1076 (2006). 

25 A. Ohtomo, D. A. Muller, J. L. Grazul, and H. Y. Hwang, 
Nature 419, 378 (2002). 

26 H. Y. Hwang, Y. Iwasa, M. Kawasaki, B. Keimer, N. Na- 
gaosa, and Y. Tokura, Nature materials 11, 103 (2012). 

27 A. Roy and D. Vanderbilt, Physical Review B 83, 134116 

( 2011 ). 

28 S. Z. Li, C. Zhu, Z. B. Yan, S. J. Luo, K. F. Wang, and 
J.-M. Liu, Journal of physics. Condensed matter : an In¬ 
stitute of Physics journal 22 , 206005 (2010). 

29 H. Chen, A. J. Millis, and C. A. Marianetti, Physical 
Review Letters 111, 116403 (2013). 

30 A. Muoz, J. A. Alonso, M. T. Casais, M. J. Martnez-Lope, 
and M. T. Fernndez-Daz, Journal of Physics: Condensed 
Matter 14, 8817 (2002). 

31 A. Azad, S.-G. Eriksson, A. Mellerga rd, S. Ivanov, J. Erik- 
sen, and H. Rundlof, Materials Research Bulletin 37, 1797 
( 2002 ). 

32 D. Iwanaga, Y. Inaguma, and M. Itoh, Materials Research 
Bulletin 35, 449 (2000). 

33 C. Khattak, D. Cox, and F. Wang, Journal of Solid State 
Chemistry 17, 323 (1976). 

34 W. Zhong and D. Vanderbilt, Phys. Rev. Lett. 74, 2587 
(1995). 

35 A. M. Abakumov, R. Erni, A. A. Tsirlin, M. D. Rossell, 
D. Batuk, G. Nenert, and G. V. Tendeloo, Chemistry of 




8 


Materials 25, 2670 (2013). 

36 R. D. Shannon, Acta Crystallographica Section A 32, 751 
(1976). 

37 G. Kresse and J. Hafner, Physical Review B 47, 558 (1993). 

38 G. Kresse and J. Hafner, Physical Review B 49, 14251 
(1994). 

39 G. Kresse and J. Furthmiiller, Computational Materials 
Science 6, 15 (1996). 

40 G. Kresse, Physical Review B 54, 11169 (1996). 

41 P. E. Blochl, Physical Review B 50, 17953 (1994). 

42 G. Kresse, Physical Review B 59, 1758 (1999). 

43 V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, 
Journal of Physics: Condensed Matter 9, 767 (1997). 

44 H. Park, A. J. Millis, and C. A. Marianetti, 

arXiv: 1505.00227 (2015). 

45 H. J. Monkhorst and J. D. Pack, Physical Review B 13, 
5188 (1976). 

46 P. E. Blochl, O. Jepsen, and O. K. Andersen, Physical 
Review B 49, 16223 (1994). 

47 A. Togo, F. Oba, and I. Tanaka, Physical Review B 78, 
134106 (2008). 

48 R. D. King-Smith and D. Vanderbilt, Physical Review B 
47, 1651 (1993). 

49 R. Resta, Reviews of Modern Physics 66, 899 (1994). 

50 G. Henkelman, B. P. Uberuaga, and H. Jonsson, The Jour¬ 
nal of Chemical Physics 113, 9901 (2000). 

51 K. Momma and F. Izumi, Journal of Applied Crystallog¬ 
raphy 41, 653 (2008). 

52 J. B. Neaton, C. Ederer, U. V. Waghmare, N. A. Spaldin, 
and K. M. Rabe, Physical Review B 71, 014113 (2005). 

53 R. Xu, S. Liu, I. Grinberg, J. Karthik, A. R. Damodaran, 
A. M. Rappe, and L. W. Martin, Nature materials 14, 79 
(2015). 

54 P. Ravindran, R. Vidya, A. Kjekshus, H. Fjellva g, and 
O. Eriksson, Physical Review B 74, 224412 (2006). 


55 P. Ghosez, J.-P. Michenaud, and X. Gonze, Physical Re¬ 
view B 58, 6224 (1998). 

56 A. Onsten, M. Ma nsson, T. Claesson, T. Muro, T. Mat¬ 
sushita, T. Nakamura, T. Kinoshita, U. O. Karlsson, and 
O. Tjernberg, Physical Review B 76, 115127 (2007). 

57 V. I. Anisimov, I. V. Solovyev, M. A. Korotin, M. T. 
Czyzyk, and G. A. Sawatzky, Phys. Rev. B 48, 16929 
(1993). 

58 M. T. Czyzyk and G. A. Sawatzky, Phys. Rev. B 49, 14211 
(1994). 

59 H. Park, A. J. Millis, and C. A. Marianetti, Physical Re¬ 
view Letters 109 , 156402 (2012). 

60 D. I. Bile and D. J. Singh, Physical Review Letters 96, 
147602 (2006). 

61 K. Butler, J. M. Frost, and A. Walsh, Energy Environ. 
Sci. 8, 838 (2014). 

62 E. Manousakis, Physical Review B 82, 125109 (2010). 

63 M. Eckstein and P. Werner, Physical Review Letters 113 , 
076405 (2014). 

64 R. Guo, L. You, Y. Zhou, Z. S. Lim, X. Zou, L. Chen, 
R. Ramesh, and J. Wang, Nature communications 4, 1990 
(2013). 

65 S. Y. Wu, W. J. Takei, and M. H. Francombe, Ferro- 
electrics (2011). 

66 Y. S. Oh, X. Luo, F.-T. Huang, Y. Wang, and S.-W. 
Cheong, Nature materials 14 , 407 (2015). 

67 J. M. Frost, K. T. Butler, F. Brivio, C. H. Hendon, M. van 
Schilfgaarde, and A. Walsh, Nano letters 14 , 2584 (2014). 

68 F. Zheng, H. Takenaka, F. Wang, N. Z. Koocher, and 
A. M. Rappe, The Journal of Physical Chemistry Letters 
6, 31 (2015). 

69 M. A. Green, A. Ho-Baillie, and H. J. Snaith, Nature 
Photonics 8, 506 (2014). 



arXiv:1506.08855vl [cond-mat.mtrl-sci] 29Jun2015 


Supplementary Information for 

’New class of planar ferroelectric Mott insulators via first principles design’ 

Chanul Kim, Hyowon Park, and Chris A. Marianetti 


1. GROUND STATE STRUCTURES 
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FIG. 1: | Fully relaxed structure of AABBOe (a) PbLaCuVOe (b) BaBiNiVe (c) BaLaCuVOe ; The colors of 
chemicla formula and structure correspond the species of atoms 
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2. ATOMIC COORDINATES, LATTICE PARAMETERS and BORN EFFECTIVE CHARGES 


TABLE I: Born Effective Charges and reduced coordinates of BaBiCuVOe in the orthorhombic 20-atoms cell. 






a = 5.387 

b = 5.457 

o 

II 

0° 

00 

o 

Atom 

7* 

rj * 

n yy 

ZL 

X 

y 

z 

Ba 

2.834 

2.809 

3.564 

0.25514276 

0.24534037 

0.26140472 

Ba 

2.834 

2.809 

3.564 

0.74485741 

0.74534082 

0.26140470 

Bi 

5.054 

4.961 

3.888 

0.26351321 

0.19642385 

0.70713054 

Bi 

5.054 

4.962 

3.888 

0.73648696 

0.69642405 

0.70713040 

V 

6.518 

4.661 

5.338 

0.24516770 

0.71660716 

-0.00331922 

V 

6.509 

4.664 

5.344 

0.75483237 

0.21660449 

-0.00331904 

C 

2.089 

2.151 

1.791 

0.24564854 

0.73988309 

0.50548329 

Cu 

2.089 

2.151 

1.791 

0.75435144 

0.23988311 

0.50548335 

0 

-3.983 

-2.191 

-1.333 

0.02899101 

0.05405394 

-0.02243997 

0 

-3.983 

-2.191 

-1.333 

-0.02899188 

0.55405366 

-0.02243998 

0 

-3.203 

-2.698 

-1.18 

0.53612190 

-0.01001588 

0.00951961 

0 

-3.203 

-2.698 

-1.18 

0.46387841 

0.48998308 

0.00951949 

0 

-1.830 

-1.913 

-3.98 

0.22966860 

0.77063504 

0.21373885 

0 

-1.830 

-1.912 

-3.98 

0.77033186 

0.27063572 

0.21373930 

0 

-2.573 

-2.584 

-2.305 

0.00094419 

-0.00745141 

0.53173795 

0 

-2.572 

-2.586 

-2.302 

-0.00094471 

0.49254838 

0.53173805 

0 

-2.562 

-2.546 

-2.502 

0.49446762 

0.00155895 

0.51652627 

0 

-2.563 

-2.547 

-2.503 

0.50553274 

0.50155854 

0.51652629 

0 

-2.347 

-2.655 

-3.275 

0.29830525 

0.79296632 

0.78021798 

0 

-2.347 

-2.654 

-3.279 

0.70169458 

0.29296668 

0.78021739 


TABLE II: Born Effective Charges and reduced coordinates of PbLaCuVOe in the orthorhombic 20-atoms cell. 






a = 5.387 

b - 5.309 

c= 8.053 

Atom 

7* 

^xx 

a 

* 

N 


X 

y 

z 

Pb 

3.575 

3.591 

4.189 

0.22158575 

0.24809749 

0.28782581 

Pb 

3.573 

3.592 

4.2 

0.72158674 

0.75190224 

0.28782594 

La 

4.343 

4.246 

5.033 

0.23719897 

0.25386950 

0.72602372 

La 

4.34 

4.247 

5.031 

0.73721805 

0.74613092 

0.72602414 

V 

4.804 

6.169 

5.659 

0.21381777 

0.74310300 

-0.00020160 

V 

4.803 

6.167 

5.655 

0.71386053 

0.25689634 

-0.00023513 

Cu 

2.122 

2.035 

1.571 

0.24320996 

0.74998738 

0.49966060 

Cu 

2.121 

2.034 

1.571 

0.74321866 

0.25001467 

0.49966488 

O 

-3.074 

-3.14 

-1.511 

-0.03795712 

-0.05741095 

-0.03023650 

O 

-2.14 

-3.751 

-1.419 

0.07014963 

0.45086863 

0.00336712 

O 

-3.075 

-3.14 

-1.51 

0.46205099 

0.05742094 

-0.03023610 

O 

-2.14 

-3.753 

-1.419 

0.57015360 

0.54911884 

0.00336638 

0 

-2.235 

-1.962 

-3.785 

0.27643088 

0.78697212 

0.21341166 

0 

-2.234 

-1.961 

-3.788 

0.77641671 

0.21303226 

0.21340697 

0 

-2.481 

-2.657 

-2.98 

0.01499215 

0.01902553 

0.51728328 

0 

-2.598 

-2.309 

-3.056 

-0.02317921 

0.51912615 

0.50752395 

0 

-2.48 

-2.657 

-2.979 

0.51498266 

-0.01903757 

0.51728362 

0 

-2.598 

-2.309 

-3.056 

0.47680972 

0.48088314 

0.50752341 

0 

-2.305 

-2.215 

-3.696 

0.28373604 

0.72132392 

0.77533943 

0 

-2.304 

-2.215 

-3.700 

0.78371752 

0.27867545 

0.77537846 
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TABLE III: Born Effective Charges and reduced coordinates of BaBiNiVOe in the orthorhombic 20-atoms cell. 






a = 5.583 

b m 5.550 

c= 7.727 

Atom 

7* 

7* 

^yy 

z L 

X 

V 

z 

Ba 

2.839 

2.912 

3.441 

0.24717643 

0.24993781 

0.25911660 

Ba 

2.839 

2.912 

3.442 

0.74717600 

0.75006205 

0.25911615 

Bi 

5.107 

4.93 

3.626 

0.19372302 

0.25087924 

0.70519973 

Bi 

5.107 

4.93 

3.627 

0.69372264 

0.74912477 

0.70519960 

V 

4.876 

6.382 

5.903 

0.21414132 

0.75021351 

-0.00276106 

V 

4.875 

6.381 

5.903 

0.71413691 

0.24978739 

-0.00276153 

Ni 

1.836 

1.798 

2.312 

0.24720605 

0.75028751 

0.50068445 

Ni 

1.835 

1.798 

2.312 

0.74720523 

0.24971204 

0.50068463 

0 

-2.551 

-3.609 

-1.265 

0.01634107 

-0.00674227 

-0.00141176 

0 

-2.502 

-3.663 

-1.289 

0.02058760 

0.50323374 

-0.00626253 

0 

-2.551 

-3.609 

-1.265 

0.51633993 

0.00674113 

-0.00140971 

0 

-2.502 

-3.663 

-1.289 

0.52058668 

0.49676624 

-0.00626069 

0 

-1.824 

-1.729 

-4.191 

0.26490638 

0.74696919 

0.22990260 

0 

-1.822 

-1.729 

-4.191 

0.76491031 

0.25303069 

0.22990495 

0 

-2.568 

-2.593 

-2.414 

0.00826344 

0.00668110 

0.52036953 

0 

-2.666 

-2.608 

-2.356 

-0.00031312 

0.50182921 

0.52601810 

0 

-2.569 

-2.593 

-2.413 

0.50826473 

-0.00667956 

0.52036855 

0 

-2.665 

-2.608 

-2.356 

0.49968819 

0.49816922 

0.52601719 

0 

-2.549 

-1.822 

-3.762 

0.28796856 

0.75605814 

0.76914401 

0 

-2.551 

-1.821 

-3.768 

0.78796864 

0.24393882 

0.76914124 


TABLE IV: Born Effective Charges and reduced coordinates of BaLaCuVOe in the orthorhombic 20-atoms cell. 






a = 5.426 

b = 5.364 

C— 8.225 

Atom 

7* 

^xx 

7* 

^yy 


X 

V 

z 

Ba 

2.704 

2.712 

3.625 

0.25296977 

0.24241584 

-0.02357038 

Ba 

2.705 

2.712 

3.628 

0.75294014 

0.74421133 

-0.02354220 

La 

4.34 

4.272 

4.067 

0.26138195 

0.24397095 

0.53009295 

La 

4.341 

4.27 

4.066 

0.76141447 

0.74299694 

0.53008170 

V 

4.486 

5.887 

5.126 

0.29296837 

0.75169868 

0.25690697 

V 

4.481 

5.887 

5.121 

0.79300514 

0.23483067 

0.25685807 

Cu 

2.149 

2.06 

1.491 

0.25227739 

0.74328231 

0.73810075 

Cu 

2.148 

2.06 

1.496 

0.75223399 

0.24358012 

0.73811651 

O 

-1.764 

-1.725 

-3.872 

0.23399499 

0.73932052 

0.04451123 

O 

-1.761 

-1.722 

-3.871 

0.73376116 

0.24739132 

0.04448562 

O 

-2.057 

-3.496 

-1.238 

-0.06691756 

-0.05553194 

0.26510464 

O 

-2.872 

-2.937 

-1.227 

0.04909435 

0.42648228 

0.26581048 

O 

-2.06 

-3.496 

-1.237 

0.43313158 

0.04201058 

0.26529792 

O 

-2.87 

-2.938 

-1.227 

0.54901425 

0.55995916 

0.26606564 

0 

-2.125 

-2.06 

-3.117 

0.21916248 

0.73796316 

0.47206989 

0 

-2.121 

-2.06 

-3.12 

0.71969351 

0.24841336 

0.47209169 

0 

-2.418 

-2.276 

-2.41 

0.00758795 

0.00039971 

0.72461584 

0 

-2.414 

-2.418 

-2.429 

-0.00266743 

0.49627268 

0.72621931 

0 

-2.419 

-2.275 

-2.408 

0.50755743 

-0.01351650 

0.72457668 

0 

-2.416 

-2.418 

-2.425 

0.49739609 

0.49064891 

0.72610664 
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3. DENSITY OF STATES 


10 

8 

6 

4 

2 

0 

10 

8 

6 

4 

2 

0 

10 

8 

6 

4 

2 


w (eV) 

FIG. 2: | Atomic resolved density of states of AABB'Oe in the R^ =7rj7rj0 structure for (a) PbLaCuVOe (b) BaBiNiVOe 
(c) BaLaCuVC>6 ; DFT + U method is employed ( Uy = 0, U cu — 5 eV, Ujvi = 5 eV); Fermi energy is set to be 0; 


(a) PbLaCuVO 


















































